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Abstract 

First-order Hamiltonian operators of hydrodynamic type were introduced by Drubrovin and 
Novikov in 1983. In 2D, they are generated by a pair of contravariant metrics g, g and a pair of 
differential-geometric objects b, b. If the determinant of the pencil g + Xg vanishes for all A, the 
operator is called degenerate. In this paper we provide a complete classification of degenerate 
two- and three-component Hamiltonian operators. Moreover, we study the integrability, by the 
method of hydrodynamic reductions, of 2+1 Hamiltonian systems arising from the structures we 
classified. 
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1 Introduction 


The theory of first-order Hamiltonian operator of differential-geometric type has been developed 
in the last three decades by several authors, starting from the pioneering work of Dubrovin and 
Novikov 1|4). In one-dimensional case, these structures are given by 


P ij =9 ij (n)^- + b^(u)u k x , 


( 1 ) 


where u = (it 1 ,..., it") are local coordinates depending on x, i, j, k = 1,... ,n, and u% = 
Dubrovin and Novikov proved that in the non-degenerate situation, namely det (<A) ^ 0, £[} de¬ 


fines a Poisson bracket through 


{F,G} = 


W pij 6G 

Su l Sui ’ 


if and only if g l i is a flat pseudo-Riemartnian metric and the coefficients r j k = !h'" h r are the 
Levi-Civita connection of the metric (where g irn g m j = ). Thus, in flat coordinates, any non¬ 

degenerate Hamiltonian operator lU} assumes constant form. In the case where the metric g is 
degenerate, that is, det(g y ) = (), this result does not hold. Grinberg | 12| and later Bogoyavlenskij 
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HI 0 firstly investigated this case, and recently we provided a complete list of two- and three- 
component Poisson structures with degenerate metric [16|. 

First-order Hamiltonian operators of differential-geometric type naturally arise in the study of 
quasilinear systems (systems of hydrodynamic type). In 1+1 dimensions they are given by 

u\ + Vj(u)ul = 0. 

Such systems are called Hamiltonian if they can be written in the form 

u\ + P ij 6jh = 0, 

where Sj = 6 / 5vP is the variational derivative, h = h( u) is the Hamiltonian density, and I nj is a 
Hamiltonian operator of hydrodynamic type 0. It was conjectured by Novikov that a combination 
of the Hamiltonian property with the diagonalizability of the matrix l/' implies the integrability. 
This conjecture was proved by Tsarev in l20l , who established the linearizability of such systems 
by the generalised hodograph transform. 

A generalisation of hydrodynamic type systems in 2+1 dimensions is given by 

u t + ^(ujuz + -B(u)u y = 0, (2) 

where u = (u 1 ,..., u n ), u l = u l (t, x, y) and A, B are n x n matrices. Systems of this type describe 
many physical phenomena. In particular, important examples occur in gas dynamics, shallow 
water theory, combustion theory, nonlinear elasticity, magneto-fluid dynamics, etc. A system 10 is 
called Hamiltonian if it can be written in the form u\ + P IJ d :i h = 0, where P IJ is a 2D first-order 
Hamiltonian operator of differential-geometric type, namely 

^ = 9lJ (U) ^ + + ~ 9lJ(u) Ty + (u K- (3) 

In 2+1 dimensions, a quasilinear system is said to be integrable if it can be decoupled in in¬ 
finitely many ways into a pair of compatible m-component one-dimensional systems in Riemann 
invariants [6j. Ferapontov and Khusnutdinova proved that the requirement of the existence of suf¬ 
ficiently many m-component reductions provides an effective classification criterion. This method 
of hydrodynamic reductions, which is a natural analogue of the generalised hodograph transform 
in higher dimensions, leads to finite-dimensional moduli spaces of integrable Hamiltonians. 

The purpose of this paper is two-fold. Starting from the classification of degenerate brackets in 
ID, we want to describe degenerate Hamiltonian operators of hydrodynamic type in 2D, that is, 
operators of the form 0 such that det(g + A g) = 0 holds V A (precise definition follows). Our analy¬ 
sis leads to a complete classification of two- and three-component degenerate structures (Section0. 
Secondly, we study the integrability, by the method of hydrodynamic reductions, of Hamiltonian 
systems arising from three-component structures we classified (Section0. 

2 Degenerate Hamiltonian operators in 2D 

The problem of classification of multidimensional Hamiltonian operators was proposed by Dubrovin 
and Novikov in J5), and thoroughly investigated by Mokhov jfT3, [l_4 ]. Some results in the classifi¬ 
cation of 2D non-degenerate Hamiltonian operators were recently obtained in our paper |9j- 

A first-order multidimensional Hamiltonian operator of differential-gemetric type (Dubrovin- 
Novikov type) is defined by 

P P = + $>)«£<*> (4) 
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where u = (u 1 ,... ,u n ) are local coordinates on a certain smooth n -dimensional manifold M or a 
domain of M n , and x = (a; 1 , x d ) are independent variables. 

As in the one-dimensional case, the condition of skew-symmetry and the Jacobi identity for a 
Hamiltonian operator of the form ((4]) impose very severe restrictions on the coefficients g' lqa (u) and 
b z f a (u). In particular, Mokhov proved the following general statement: 

Theorem 1 ( lfl5ln . Operator of the form ((Tj) is a Hamiltonian operator, i.e. it is skew-symmetric and satisfies 
the Jacobi identity, if and only if the following relations for the coefficients of the operator arefidfilled: 


S y “ = 9 3ia , 


dg ija 

du k 


=&r+&r, 


E 

, dbj r 


J2 ( g sia V s r > 3 - g sj0 bi ra ) = 0, 

(«./?) 

Y {g sia K p ~ g sjp bf a ) = o, 

13b> r P db jrp \ 

„sia I s _ g | _|_ ^ija^srp _ ^irot^sjp 


du q du s 


dbj , 1 


« ”q 


= 0, 


si/3 —- b¥ p bi ra - bi r Pbj sa = g sja ^f - b>f a b 8rp - 6“^K r “, 

y du s 8 q 8 q y du s 8 q 98 


(5a) 

(5b) 

(5c) 

(5d) 

(5e) 

(5f) 


d 

du k 


' BK> rp db jr0 \ 

UUs _ q I i hijot^srP _ virai sjP 

du q du s 8 9 8 9 


8 

du q 


si/3 


dbj r 


du k 


db{ ra 

du s 


b ijp b s k ra - bl rl3 b 8 k J 



= 0. (5g) 


Relations lf5al and d5bt are equivalent to the skew-symmetry of the bivector 0, and relations 
(l5cl— (5g i are equivalent to the fulfilment of the Jacobi identity for a skew-symmetric bivector of the 
form 0. The signs ^ and ,. . fc . mean cyclic summation on the indicated indices. Notice 
that for n = 1, these conditions reduce to Grinberg's conditions |l2j. 

In the one-dimensional case, Hamiltoinan operator 0 is called degenerate if det (g iq ) = o. In the 
multidimensional situation, we have the following 


Definition 1. A d-dimensional operator of the form 0 is said to be degenerate if det (Ea=i = 0 

for any choice of X\,... ,Xd, be. if there is no linear combination of the metrics g a such that the determinant 
of this linear combination is non-zero. 


Let us point out that Theorem [T| does not assume non-degeneracy of operators or additional 
conditions on the coefficients of 0. 

From Mokhov's conditions it immediately follows that each multidimensional Hamiltonian op¬ 
erator of the form 0 is always the sum of one-dimensional Hamiltonian operators with respect to 
each of the independent variables x a , fl5| . 

Based on this result, and on the classification of one-dimensional degenerate Poisson structures 
of hydrodynamic type, we give a complete description of two- and three-component degenerate 
Hamiltonian operators for d = 2, namely 

P ij = g l3 ( u)E + b£(u)u k x + u)E + b£{u)u k . (6) 
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For simplicity, let us label the cc-part and the y -part of the Hamiltonian operator © respectively 
with P( x) and P (j/) . 

Hamiltonian structure of the form 0 is called trivial if it is identically zero, or if it can be reduced 
to the form 

~9 ij =t9 ij , %=&*, (7) 

for £ constant. Notice that allowing linear change of the independent variables x, y, an operator 
satisfying 0 is essentially ID. 

Remark. Let us remark that if a pair of Hamiltonian operators defines a 2D structure, by 0 it 
easily follows that these two operators are compatible, and therefore they define a bi-Hamiltonian 
structure fT4lll5ll . Degenerate bi-Hamiltonian structures of hydrodynamic type were firstly inves¬ 
tigated by Strachan BT8l[T9fl , revealing a nice relation with the theory of the analogous of Frobenius 
manifolds with degenerate metric. 


2.1 Classification 


The analysis of Mokhov's conditions 0 is not straightforward. In order to study two- and three- 
component structures, we fix the pair ( g , b) given by the classification of ID Hamiltonian operators 
[16). This classification can be summarised in the following two theorems. 

Theorem 2. Any degenerate two-component Hamiltonian operator ofDubrovin-Novikov type in ID can be 
brought, by a change of the dependent variables, to one of the follozving two canonical forms: 


P = 


0 


f d 

U> X i 


p = 


2 u 
U x n 

-f 0 

u L 


( 8 ) 


Theorem 3. Any degenerate three-component Hamiltonian operator of Dnbrovin-Novikov type in ID can 
be brought, by a change of the dependent variables, to one of the follozving canonical forms: 


rank(g) = 0: 


P = 



(9) 


rank(g) = 1: 



^ d,- 

0 


^ dx 

*4 

°) 

^ d x 

0 

_«2\ 
u 1 \ 

^ d x 

u 1 

-4\ 

U 1 \ 

p = 

0 

0 

o\,P = 

~ul 

0 

0 , P = 

0 

0 

o 

II 

U 2 X 

u^- 

0 

0 



0 

«) 

K 0 

0 

oj 

u 

\ U L 

0 

0 ) 

14 

N U L 

0 

0 ) 


( 10 ) 


rank(g)=2: 


P = 


/o 

dx 

°\ 


( 0 

dx 


0 

0 

,p= 

dx 

0 

v° 

0 

oj 


1 < 
\7A 

0 


,P = 


0 

,, 3„,3 


U°U^ — U‘ 


U°U X —U* 


0 


( 11 ) 


Once we have fixed the pair (g, b), solving 0 we are able to find the pair (Jj. b). At this point, 
we look for canonical forms of 2D structures using transformations which preserve the form of the 
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first structure given by (g, b). As we will see, in some cases these transformations are not enough 
to eliminate all the functional parameters appearing in the 2D structure. 

Let us agree on some notation: if a function depends only on one variable, we denote with ' 
the derivative with respect to that variable. Otherwise, if a function depends on more than one 
variable, say, / = /(it 1 ,..., u n ), then we use dif = jjj-. In Section[3j for simplicity, the derivative 
with respect to id will be denote as f u <. 


2.1.1 Two-component case 

Here we provide a full description of the two-component case. 

Theorem 4. Any non-trivial degenerate two-component Hamiltonian operator of Dubrovin-Novikov type 
in 2D can be brought, by a change of the dependent variables, to the following form 


u _ (d x + u 2 d y + \u 2 -e 
± — I 2,22 

I U +U M 

\ e—m— ~ 


<+u'K 


( 12 ) 


where e can be either 0 or 1. 


Proof: 


First of all, the case g = g = 0 gives no non-trivial solutions. In the case where the rank of the pencil 
</® J ’ + Xg v is constantly equal to one, there exists a coordinates system (u 1 , it 2 ) where 



/ O' 

0 o, 


here / = /(it 1 , u 2 ) is some function. Let us fix the //.,.) structure. 

Case (8) 1 . If K are all identically zero, conditions (0 imply 

f = f(u 2 ), = y, 

and all other b l f equal to zero. If / = £ is constant, than g = £g. Otherwise, using a transformation 
which preserves P( x ), that is, a suitable change of coordinates of the form u 1 = v 1 + ip 1 {v 2 ), u 2 = 
ip 2 (v 2 ), we can easily reduce / to v 2 , obtaining lfl2ll with e = 0. 

Case {8} 2 . Suppose b 21 = —bf 2 = dp. Conditions (0 imply 

/ = /(« 2 ). ^ = t? = -i.J 2 = L. 

If / = ^ is constant, than g = £g, b = £b. Otherwise, let us assume / non-constant. Transformations 
which preserve are given by u 1 = v 1 , u 2 = y>(v 2 ), then we can always choose such that / 
reduces to v 2 in the new coordinate system, obtaining dl2l) with e = 1. ■ 


2.1.2 Three-component case 

The analysis of the three-component situation is more complicated. Let us consider separately the 
cases with respect to the rank of the pencil g\=g — A g. We point out that in some cases the group 
of transformations preserving the first structure \\ x) is not sufficient to reduce the second structure 
Ify , to something simpler. Thus, we will just consider the more general structure given by the 
solution of Mokhov's conditions. The results can be stated as follows. 
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Theorem 5. Rank(g>) = 0. Any non-trivial degenerate three-component Hamiltonian operator of Dubrovin- 
Novikov type in 2D can be brought, by a change of the dependent variables, to one of the following forms: 



/ 

0 

ul+ifv 3 


/ 

0 

u 3 x + U 3 U 3 

°\ 

p = 

~U 3 x 

— U l Uy 

0 

o 

73 

II 

-4 

- u 3 ul 

0 

° 


V 

0 

0 

«/ 

\ 

0 

0 



(13) 


In this case, we do not need any linear change of the independent variables x, y. 


Theorem 6. Rank(g>) = 1. Any non-trivial degenerate three-component Hamiltonian operator of Dubrovin- 
Novikov type in 2D can be brought, by a change of the dependent variables and linear change of x and y, to 
one of the following forms: 


( d x 


P = 


(U 2 dy A 


— hlly 


hu v 

0 

0 


( d x 


,P = 


fdy 


fofu 2 +d 3 fu 3 


0 


u^—hvf+fvC: 


u l~ hu l+f u l \ 

- TP - 

0 

0 ) 


(14) 



(d x +id v +^<* s ’>< 

u 3 x + hu 3 

°\ 

( d x 

+ Vfdy + \ 

2 . 2 2 
u x +u u y 

u 1 

u 3 +ll 2 

u 1 


p = 

-v 3 - hu 3 y 

0 

0 ,P = 


u 2 x + u 2 u 2 y 

U 1 

0 

0 



\ o 

0 

0/ 

V 

u 3 x +u 2 ul 
u 1 

0 

0 

/ 


(15) 


where f = f{u 2 , u 3 ), h = h(u 2 , u 3 ) are arbitrary functions and e can be either 0 or 1. 


Theorem 7. Rank(g>) = 2. Any non-trivial degenerate three-component Hamiltonian operator of Dubrovin- 
Novikov type in 2D can be brought, by a change of the dependent variables and linear change of x and y, to 
one of the following forms: 


( —2 u 1 d v — u}. d x + u 2 d v + 2 u 2 


P = 


*y 

9 

' U 

—eu 


y 

y -u: 

3 

y 


y 


( 


p = 


P dy 


r_% 

2 


dx + q dy + (q' - e)i 


V 


0 


d x + q d y + eu 3 
r d y + —^ L 
0 



P = 


P = 



P = 


ed 


•y 

T U 3 dy | V,y 
ul+u 3 u 3 v 

- TV —" 


■ U7 


dx + U 3 dy — 

0 

0 


u x +u 3 ul \ 


P = 


d x 




P = 


0 

d x 


±<\ 


\dy 


0 

0 

p = 

p = 


p = 


^ dr ~^ d„ - 


d - — d + 

u* x 2 ' 2 

,,3 


V 


2 

0 

0 


J 

-4\ 

0 

0 


dy 

d x — u 3 d y — u 3 

\ u 3 x -2 u 3 u 3 y 

\ u^u 1 —u 2 

/ K, dy 


d x 


ut- 2u* 


U 


Xy 


Ajy I il ILy 


(u 3 ) 2 dy + U 3 U 3 — 


TTTfTr 


u—u* 


KU y 

T 77T\ ! I 


d X ~Kdy 
\ U x~^ KU y 

\ u'^u 1 —u 2 


dx — Kdy + 

KV? cL + - 


u ut —2 ku ut 


u^u 1 —u 2 


u 3 ul- 2(u 3 ) 2 

< 

vPu 1 —u 2 


0 

J 

ul-2K.nl 

\ 

u°u x —u* 

u 3 u 3 - 2 ku 3 ut 

—rrm—— 


0 


J 


where p, q, r are arbitrary functions on u 3 , k is constant and e can be either 0 or 1. 


(16) 

(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 
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The proof of these theorems can be found in the Appendix. 

Let us point out that, after swapping the coordinates u 1 , u 2 , dTSb o corresponds to the Hamilto¬ 
nian operator for the 2D equations of gas dynamic (see, for instance, f7j), namely 


pij _ 


(o 

dx 

\dy 


dy A 


u 2 -ul 
y x 


We will discuss it in Section lT2l 


( 22 ) 


3 Hamiltonian systems of hydrodynamic type in 2+1 dimensions 

In this section we discuss (2+l)-dimensional Hamiltonian systems of hydrodynamic type, 

u t + A(u)ua; + -B(u)u.y = 0, (23) 

which are representable in the form u t + Ph u = 0, where h{ u) is a Hamiltonian density and P 
is a two-dimensional Hamiltonian operator of differential-geometric type l(6|. As we recalled in 
the introduction, a (2+l)-dimensional quasilinear system is said to be integrable if it can be decou¬ 
pled in infinitely many ways into a pair of compatible m-component one-dimensional systems in 
Riemann invariants. Let us briefly describe the method of hydrodynamic reduction introduced by 
Ferapontov and Khusnutdinova in [6], 


3.1 The method of hydrodynamic reductions 

The method of hydrodynamic reductions is based on the existence of exact solutions of the (2+1)- 
dimensional system (f23l> . These solutions have the form u = will 1 ...., R m ), where the Riemann 
invariants R = (R 1 ,..., R m ) solve a pair of commuting diagonal systems 

R\ = A ,; (R) /?;, P! y = m‘(R) i?;. (24) 


Let us point out that we do not impose any constraint on the number of Riemann invariants: m 
is arbitrary. Therefore, the (2+l)-dimensional system we are considering (f23l> . is decoupled into a 
pair of diagonal (l+l)-dimensional systems given by ll24t . Usually, these solutions are known as 
nonlinear interactions of m planar simple waves. 

It turns out that the commutativity of the flows d24t is equivalent to the following constraints 
on the characteristic speeds A*, (20l : 


_ djV g - d 

X> - A i fp - T 3 dRi ’ 


(25) 


(no summation). Imposing these restrictions, the general solution of systems (l24ll is given by the 
implicit generalised hodograph formula [20] 


t/(R) = x + A*(R) t + R) y, i = 1,..., m. 


(26) 


Here the functions v l (R) are characteristic speeds of the general flow commuting with (l24b . namely, 
the general solution of the linear system 


dnV <9,A dj/j, . 

- — -7 = - — -7 = - — -7, U). 

vl - v l \i - A* li? - 


(27) 


7 












By straightforward computation, the substitution of u (R 1 ,..., R m ) into (123b , using < f24b . leads to 

(E + AM + fj, l B)diU = 0, i = (28) 

where E is the n x n identity matrix. This means that both A* and //' have to satisfy the dispersion 
relation 

det(.E + X i A + fSB) = 0. (29) 

Furthermore, the construction of nonlinear interactions of m planar simple waves can be sum¬ 
marised as follows. First of all, we have to decoupled the initial (2+l)-dimensional system J23b into 
a pair of commuting flows d24b , by solving the equations (25) , (28) for u(R), A ! (R), M’(R) as func¬ 
tions depending on the Riemann invariants R 1 ,R m . It is not difficult to see that for m > 3 the 
system given by these equations is overdetermined. Thus, in general this system does not posses 
solutions. However, if we are able to construct a particular reduction of the form (24), the sec¬ 
ond step is quite straightforward: we have to solve the linear system given by (27) for the functions 
■(/' (R), and then we can obtain R 1 ,R m as functions off. x. y from the implicit hodograph formula 
(26). 

What can we say about the number of //(-component reductions that a (2+l)-dimensional system 
(23) may admit? Analysing equations (25) and (28) . one can prove that this number is parametrised, 
up to changes of variables of the form R 1 —>• /’(/?' ), by to arbitrary functions of a single variable. 
Remarkably, this number does not depend on n. This leads to the following definition. 

Definition 2 (03). A (2+l)-dimensional quasilinear system is said to be integrable if it possesses m- 
component reductions of the form (24) parametrised by m arbitrary functions of a single argument. 

Remark. Looking at the structure of equations (25) and (28) , one can see that their consistency 
conditions involve only triple of indices i ^ j k. Moreover, all these conditions are completely 
symmetric ini, j and k, and then it is enough to verify them setting, for instance, i = 1, j = 2, k = 3. 
This means that the existence of non-trivial three-component reductions implies the existence of m- 
component reductions for arbitrary to. 

Remark. We require that A* and p‘ do not satisfy any linear relation, otherwise we would have no 
sufficiently many arbitrary functions of a single argument. Indeed, let us suppose that p l = a A' + b. 
Condition (25) reads djaX 1 + djb = 0, which implies a and b constant. Thus, solutions of the system 

m, 

R\ = A‘i£, R l y = ( a\ l + b)Rl 

are of the form R 1 = R‘(x + by,t + ay). These solutions correspond to travelling wave reduction, 
and they clearly do not contain enough arbitrary functions. 

3.2 Generalised two-dimensional gas dynamic equations 

The equations of two-dimensional isentropic gas dynamics are of the form 

Px P 

Pt + ( pu) x + (pv)y =0, lit T UU X + VUy + = 0, v t + uv x + VVy + = 0, (30) 

where p = p(p) is the equation of state. In matrix form (23) , one has u = (p, u, v) f and 



/ u p 0 \ 


/ v 0 p \ 

A = 

c 2 

— u 0 

, B = 

0 v 0 
ji 


3 

o 

CL O 


c 

7 / 
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where c 2 = p'(p) is the sound speed. As demonstrated in fl7j , there exist potential flows describing 
nonlinear interaction of two sound waves which are locally parametrised by four arbitrary func¬ 
tions of a single argument. 

The system <f30l) can be written in Hamiltonian form as u* + Ph u = 0, where the operator P is 
given by ll22ll . namely 



the Hamiltonian density h is h(p, u, v ) = \p(u 2 + v 2 ) + k(p), and the equation of state p and the 
function k are related by p p = pk pp . 

Let us assume h = h(p , it, v) generic, thus the system u t + Ph u = 0 reads 

qi _ qj qj — qi 

Pt + ( h u )x + (h v )y = 0, Ut + {hp) x H—— h v = 0, v t + ( h p ) y H -— —h u = 0. (31) 

Let us consider the Riemann invariants R 1 ,..., R m solving 

Ri=X i (R)Ri R;=p\R)Rl, i = 1,... ,m. 

By straightforward computation, the substitution p = p(R), u = it(R), v = i>(R) into d3it implies 

(1 + A l h pu + p i h pv )diU + h pp Xdip = 0, (32) 

(1 + A l h pu + p l h pv )diV + h pp p l d i p = 0, (33) 

(1 + A ‘h pu + p!h pv )dip + (X i h uu + p‘h uv )diU + (A i h uv + p i h vv )diV = 0, (34) 

here i = 1,..., m, di = Note that since p l diU = A l diV (this easy follows from ( f32t and J33t , 
assuming 1 + A’ h fm + p l h pv ^ 0), one has u y = v x . Thus, solutions are necessarily potential. Then, 
setting u = ip x and v = p y , our system d3Tb reads 

Pt T ( ku)x T (kv)y = 0; <Pxt “b {kp)x — 0, Pyt “t“ ( h p )y = 0- (^^) 

Both the last two equations give pt + h p = 0, so we finally have the following system 

Pt + (h u )x + ( h v ) y = 0, ipt + h p = 0. (36) 

If we consider the partial Legendre transform 

p = h p , u = u, v = D, h = h — ph p , (37) 

the derivatives respect the new variables are 

hp - P, hu - h Ul hy - hyi (^^) 

and we can rewrite the system ( f36l> in the form 

(hp) t + (hu) x + (hy)y = °, <Pi=/5, <fix=u, ip v = v. 

The function h depends only on p x . ip y , <p t and thus we obtain three-dimensional Euler-Lagrange 
equations (setting h = f) 

(UX + (Uy)y + (Mt = o> ( 39 ) 
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corresponding to Lagrangian densities of the form f(<p x , ip y ,(pt)- For example, the Lagrangian den¬ 
sity f = u^ + Uy — 2e Ut leads to the Boyer-Finley equation u xx + u yy = e Ut u t t 0. 

In |8) Ferapontov, Khusnutdinova and Tsarev derived a system of partial differential equations 
for the Lagrangian density f((p x ,ip y ,ipt) which are necessary and sufficient for the integrability 
of the equation (f39l) by the method of hydrodynamic reductions (see also ItlOl for further details). 
Setting a = ip x ,b = (p y ,c = ipt, these conditions can be represented in a remarkable compact form: 


is integrable by the 


(40) 


Theorem 8 (US)). For a non-degenerate Lagrangian, the Euler-Lagrange equation 
method of hydrodynamic reductions if and only if the density f satisfies the relation 

d A f = d3 f-ff- + jf det(dAf); 

here d 3 f and d A f are the symmetric differentials of f. The Hessian H and the 4 x 4 matrix M are defined as 
follows: 

( 0 f a 

fa 

h 
V fc 


H = det 


^ faa fab 

fab fbb 

y fac fbc 


M = 


faa 

fab 

fac 


fb 

fab 

fbb 

fbc 


fc 

fac 

fbc 

fee 


\ 




The differential dM = M a da + M^db + M c dc is a matrix-valued form 


dM = 


/ 0 

faa 

fab 

fac ^ 


( 0 

fab 

fbb 

fbc 

\ 

faa 

faaa 

faab 

faac 

da + 

fab 

faab 

fabb 

fabc 


fab 

faab 

fabb 

fabc 

fbb 

fabb 

fbbb 

fbbe 


\ fac 

faac 

fabc 

face j 


V fbc 

fabc 

fbbe 

fbcc 

/ 



( 0 

fac 

fbc 

fee N 






+ 

fac 

faac 

fabc 

face 

dc. 






fbc 

fabc 

fbbe 

fbcc 







\ fee 

face 

fbcc 

fccc ) 






db 


(41) 


Finally, we recall that the equations of gas dynamic possess only double waves reduction, and 
are not integrable by the method of hydrodynamic reductions [7J. On the other hand, the gener¬ 
alised equations Oil define a (2+l)-dimensional integrable system when the Lagrangian density 
f((p x , 'Ey ■ Et), obtained by the Flamiltonian density hip. u. v) performing a partial Legendre trans¬ 
form 071 . satisfies the conditions given by Theorem[8] 


3.3 Three-component Hamiltonian systems with degenerate structure 

We have seen that the degenerate Hamiltonian operator fill ., leads to a class of integrable systems 
related to the Lagrangian density of the form f(tp x ,Ey,Et)- Here we are going to describe all three- 
component cases arising from our classification. 

The aim of this section is to apply the method of hydrodynamic reductions to three-component 
Hamiltonian systems given by u t + Ph u = 0, where P is a Hamiltonian structure appearing in 
Theorems [5j [6] and [7] Let us identify the Hamiltonian operators we obtained with the rank of the 
pencil g\. For instance, we call rank-zero structures the Hamiltonian operators listed in Theorem[5] 

Theorem 9. The method of hydrodynamic reductions imposes additional differential constraints under 
which equations under study reduce to known classes of systems considered before: 

• rank-zero structures lead to trivial systems 

u\ = u 2 t = Ut = o, 
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• rank-one structures lead to one dimensional system of the form 


u\ + fW-)v} x = 0 , 


u t =u t = 0 , 


• rank-two structures lead either to one dimensional system to the form 

ul + (h u - 0x=0, Uf+(h u i) x =0, zzf = 0, 
or two-component non-degenerate Hamiltonian systems 

ul + {h u i) x = 0, v% + (h u 2 ) y = 0, (42) 

ul + ( K 2 ) x = 0, zz 2 + (h u i ) x + (/i „2 ) y = 0, (43) 

zz^ -(- (2zz h u i u h u 2 — h^ x (zz /z u 2 )y = 0, u^ -(- (zz h u i )^ 4 (2zz h u 2 -t- zz h u i /z)y = 0, (44) 
p/zzs f/ze trivial equation vf = 0, or to the system 

u] + (h u 2 ) x + (h u s)y = 0, vf + {h u i)x=0, rz? + (/i„i) w = 0. (45) 


We point out that the integrability of two-component non-degenerate Hamiltonian systems d42t . 
J43l l and ( 144b . generated respectively by the Hamiltonian operators 



is completely understood, see IfTTI for further details. Furthermore, as we showed above, sys¬ 
tem d45l> reduces to the three-dimensional Euler-Lagrange equations J39b after performing a partial 
Legandre transformation of the form d37b . 

Proof of Theorem[9l 


First of all, let us remark that if u\ = 0, for some i, the method of hydrodynamic reductions 
necessarily implies u 1 = const. Secondly, if one of the equations of the system is of the form 
ui+<t>(u)ui+ip(u) Uy = 0, the method of hydrodynamic reductions implies (A- 5 + f + ipp,^ )djU l = 0, 
which leads to u l = const, since we are imposing that A J and fi :i do not satisfy any linear rela¬ 
tion. Furthermore, in these cases we can replace u l with a constant, and then the Hamiltonian will 
depend on a :l for j f i. 

Using these observations, the proof is straightforward. Rank-zero structures easily lead to trivial 
systems. For the rank-one structures we always have zz 2 and zz 3 constant, which leads to an operator 
of the form 


P = 


(d x + nd y 

0 

V 0 


0 o\ 
00, 
0 0 J 


k = const , 


which is essentially one-dimensional (up to linear change of the independent variables x and y). 

The analysis of rank-two structures is a bit more complicated. In the cases | [T6t 1 and llT9L. the 
method of hydrodynamic reductions implies zz 3 = const. Thus, up to a change of local coordinates 
zz 1 , zz 2 , the 3x3 degenerate Hamiltonian operator reduces to direct sum of the 2x2 two-component 
non-degenerate Mokhov's Hamiltonian operator 143 44) 


P = 




d x + 
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and the trivial lxl operator P = 0. 

In the cases (1171 ] 2 , d l 81 , , d20b and (f21~b the method of hydrodynamic reductions implies again 
u 3 = const. These structures reduce to direct sum of constant 2x2 two-component non-degenerate 
Hamiltonian operator, and the trivial lxl operator P = 0. Constant 2x2 non-degenerate Hamil¬ 
tonian operators are known [11 j: if they do not reduce to one-dimensional operator 



(for instance, for e = 1), they can be brought to one of the following two forms 

p _ (dx 0 \ p _ ( 0 d x \ 

y 0 dyj ’ \d x dyj ’ 

by a change of local coordinates u 1 , u 2 and a linear change of the independent variables x, y. 

It remains to consider the cases (TT6b o and (1191]. It is not difficult to see that in both cases we get 
u 2 = u 3 . Therefore, solutions are necessarily potential. Then, setting u 2 = ip x and u 3 = tp y , the 
system leads to ( f35l> . ■ 

4 Concluding remarks 

The problem of classification of 2D Hamiltonian operators of differential geometric-type, proposed 
by Dubrovin and Novikov in 0, is now completely solved up to three-component case. Even 
though in f9j we provided a complete classification of non-degenerate operators up to four compo¬ 
nents, in the degenerate case it is still open. The main obstacle is the lack of a full description of one¬ 
dimensional degenerate Poisson brackets. Indeed, already for four-component one-dimensional 
degenerate structures, the computation of Jacobi conditions is quite complicated fT2lll6ll . 

As we have said above, any 2D degenerate Hamiltonian operator gives rise to a pair of ID 
compatible degenerate brackets of Dubrovin-Novikov type (see [14.151 for further details). Some 
of the degenerate bi-Hamiltonian structures arising from our classification are not of the kind in¬ 
vestigated by Strachan l H8lfl9l . It would be interesting to analyse these structures and to study a 
possible correspondence with the analogous of Frobenius manifolds with degenerate metric. 
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Appendix. Proof of Theorems HJ,[6] and [7] 

Proof of Theorem^ 

In the case where the pencil g\ has rank constantly equal to 0, both the metrics must be identically 
null. Thus, by Theorem[3j we can always reduce the coefficients bV to bl 2 = —b 21 = 1, that is, the x- 
part of the 2D Hamiltonian operator can be fixed as 0. Imposing 10, we get that all the coefficients 
must vanish except b [2 = —//]' = v(u [ . u 2 . u 3 ). The transformations which preserve the form of 
P( :r) have the form 

u 1 = <p 1 (ri 1 ,t; 2 ,'y 3 ), u 2 = ip 2 ^ 1 ,v 2 ,v 3 ), u 3 = ip 3 (v 3 ), 
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with the constraint 


diip 1 d 2 <p 2 - d 2 p l dip 2 = (<p 3 )'. 

If v = £ is constant, we get b = !pb. Let us assume that v = v(;u' > ). Then, in the new system 
of coordinates, it is always possible to reduce 1 / to w 3 . Let us finally assume that v is an arbitrary 
function of u 1 , u 2 , u 3 . Then, there exists a change of coordinate preserving P( x ) which transforms v 
to v 1 or, equivalently, to v 2 (these two cases are the same since we can swap v [ , v 2 ). Summarising, 
P( y ) leads to one of the following two structures 



/ 

0 

V 1 V 3 



/ 

0 

V 3 Vy 

<A 

P (y) = 


-v 1 v 3 

0 


> P (y) = 


-v 3 v 3 

0 

° 


V 

0 

0 

V 


V 

0 

0 



which are not equivalent modulo transformations which preserve the form of P( x )- 


Proof of TheremU] 


When the rank of the pencil g' 1 ' 3 +Xg lJ is constantly equal to one, we can always work in a coordinate 
system where the metrics assume the forms 



fl 

0 

°\ 

(f 

0 

°\ 

<f J = 

0 

0 

0 , g lJ = 

0 

0 

0 



0 

oj 


0 



where / = f(u [ . u 2 . u 3 ). Let us now consider separately each case given by Tire ore in [3] 

Case lITOlq. The symbols b'P are identically 0. In this case, a generic transformation which preserves 
P(x) is given by 

1 1 1 1 / 2 3 \ 2 2 / 2 3 \ 3 3 / 2 3 \ (A/i\ 

u =v + ip (v ,v ), u = ip {v , v ), u = ip {v , v ). (46) 

Let us point out that this change of coordinates transforms b l f! (and then If ) as components of a 
(2, l)-tensor fl6l . Conditions © imply two solutions. 


Solution 1. The first solution reads 


-f tt 2 3\ 711 9 2 f ~ n d 3 f ~21 712 713 731 / 

f = f{u,u), b 2 = —, b 3 = —, b 2 =-b 2 = bf = -b% =i/j, 


ip 2 
V 


b 3 = ~b 3 = f—, bf = -bf = 77, 


where ip = ip(u 2 ,u 3 ) and 7/ = ri(u 2 ,u 3 ), and all other fjjp vanish. Clearly, we are imposing 77 f 0. 
The operator leads to 

d2ful+d3ful 


(fd 


P (y) = 


y 


rj y 

V -r}ul - V’W 3 







where / = f(u 2 , u 3 ), ip = ip(u 2 , u 3 ) and 77 = 77 (u 2 , u 3 ) f 0. Applying a transformation of the form 
we get / -7 f(p 2 ,p 3 ) and 


-b 12 = bf = -bf = iP 


-7 



bf = -bf = 77 


31 




{d 3 p 2 rj + d 3 p 3 ip)(d 2 p 2 rj + d 2 p 3 ip) 
( d 2 p 2 d 3 p 3 - d 3 p 2 d 2 p 3 )r\ 

( d 3 p 2 r] + d 3 p 3 ip ) 2 
(d 2 p 2 d 3 ip 3 - d 3 ip 2 d 2 p 3 ) 77’ 

( d 2 p 2 rj + d 2 (p 3 ip ) 2 
(d 2 tp 2 d 3 p 3 - d 3 ip 2 d 2 ip 3 )r]' 
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We cannot choose both dip 2 r\ + ct^ 3 '*/’ = 0 and d 3 <p 1 r] + = 0, otherwise we would have 

the denominator equal to zero. However, a suitable choice of the functions p 2 and tp 3 allows us to 
reduce / to either v 2 or v 3 (which are equivalent up to swapping v 2 and v 3 ) if / is not constant, and 
if) to zero. This leads to two operators 



( nd y 

0 

VV 2 y\ 

( V 2 dy + \ 

0 

V v l\ 

•a 

II 

0 

0 

0 1 5 P (v) = 

0 

0 

° 



0 


\ -d v l 

0 

0 / 


where k is constant and ?/ = fj(v 2 , v 3 ). Allowing linear change of x and y, n can be brought to zero. 
Solution 2. In the case where 77 = 0, the solution reads 

and all other tfj 2 vanish. Modulo transformations of the form d46t this case corresponds to the 
previous one. 

Case m,. Here bl 2 = —b 21 = 1, while other symbols tfj 2 are identically 0. Conditions (0 imply 

/ = /(« 2 .» 3 ), 65 ' = ^, = iy = -i? = v(u\uy, 

and other b'J = 0. A generic transformation which preserves P( x \ is given by 

u 1 =v 1 , u 2 = d 3 tp 3 (v 3 )v 2 + p 2 (v 3 ), u 3 = ip 3 (v 3 ). 

Unfortunately, in general this group of transformations cannot help to simplify our structure (the 
operator depends on two functions of u 2 , u 3 , while the group depends only on two functions of u 3 ). 
We should consider separately each case where the functions / and v are constant or depend on 
one single variables. Therefore, it is more reasonable to consider just the general solution, namely 


P iv) = 


(fdy + d2fU V+ d3fU « UUy 0 N 
0 0 


—vu 


\ 


y 


0 


0 0 ; 


for arbitrary / = f(u 2 , it 3 ), v = v(u 2 ,u 3 ). 

Case m 3 - Here & 3 1 = ~b\ 3 = while other symbols b) 2 are identically 0. Conditions 101 imply 


j- jv 2 3 \ rn ^2/ 

f — f\ u i u )■> b 2 — ——, w 3 


b\ l = 


d 3 f 


bl 1 = -bl 3 = 


f 


bl 3 = -bV = 


;3i 

2 1 : 

U l 


where v = v(ii 2 . u 3 ), and other }>f = 0. A generic transformation which preserves I :> ( x ) is given by 

u 1 =v 1 , u 2 =ip 2 (v 2 ), u 3 = ip 3 (v 3 ). 


As before, this group of transformations cannot help to simplify our structure for arbitrary / and 
v. Therefore, the operator leads to 


/ fdy + 


P (v) = 


vu v~ 

-^T 


!A\ 
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for arbitrary / = f(u 2 ,u 3 ), v = v(u 2 ,u 3 ). 

Case dTob i. Here b 21 = —b\ 2 = — = 63 1 = p-, while other symbols b l £ are identically 0. Condi¬ 

tions © imply 


f = f(u 2 ,u 3 ), b\ l = 


82 f 


d 3 f 


‘“= 2 ' 


K' = ->>'/ = = K‘ - 


,1’ 


and other b) : ! = 0. A generic transformation which preserves P( x ) is given by 

u 1 = v 1 , u 2 = tp 2 (v 2 ,v 3 ), u 3 = ip 3 (v 2 , v 3 ). 


This change of coordinates transforms the objects b'ql as components of a (2, l)-tensor Ifl 6 l . If / = £ 
is constant, then g = £<? and b = £ 6 . Otherwise, we can choose <p 2 or <p 3 such that / reduce either to 
(r or v 3 , which are equivalent forms since we can swap v 2 , v 3 . Thus l , ( Vl leads to 


P(v) 


d + -S- — 

u, y 1 o 


*y 

2 

V V. 


0 

0 


Proof of Therem[7] 

When the rank of the pencil g l] + A g^ is constantly equal to two, we have three possibilities: 


or 



(47) 


(48) 


where p, q, r are arbitrary functions of u 1 , it 2 , u 3 . We remark that 

,,2 


and 


(49) 

are equivalent up to a 


transformation of the form u 1 = v 2 ,u 2 = v [ , v 3 = v 3 (which preserves the form of the first metric). 
However, this change of coordinate does not fix the structures CD} 2 ar *d (I I It ,. Therefore, when the 
x-part of the operator is given by (TT) ,, we can avoid J49l >, while in the other two cases, and 
m 3 , we have to take it into account. 

Case On Here the first structure is (H U and the group of transformations which preserve its 
form is given by 


u 1 = e^v 1 + ^ 5 1 (v 3 ), 


u 2 = e ^v 2 


+ y>V), u 3 = v 3 (v% 


(50) 


where d is constant, plus the switch of u 1 , u 2 (note that this change of coordinates transforms the 
objects b ^ as components of a (2, l)-tensor fl" 6 j ). In the case where we are dealing with 147k up to 
swapping u 1 , u 2 , solutions of conditions faj can be summarised as follows. 
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Solution 1. The first solution is given by 


p = p(u 3 ) + KU \ q = q(u 3 ) — b 1 1 1 = b 2 2 1 = ^, &3 1 = y, 

b 2 2 = — k, 6g 2 = 2g / , 6 21 = —q', 

where k ^ 0 is constant. This leads to 

* (p + ku 1 ) d y + KUy ^ P u v d y -KUy + 2 q'u 3 o\ 


P (y) = 


\ 


(q~¥) dy + ^-q'u: 


/„,3 

y 


0 

0 / 


A change of coordinates of the form (foOt transforms 


p -t p 1 K+p{p 3 ) 1 q 


ip 2 K 


+ q(ip 3 ), k —>■ ne ^ 


thus, it is always possible to choose y 1 , y 2 and ip such that in the new coordinates p = q = 0 and 
is fixed, let us set it equal to —2. Therefore, the operator leads to 


/ —2V 1 d y — Vy v 2 d y + 2Vy 0^ 


P (y) = 


\ 


0 


Let us point out that in this case the 2D operator P can be view as direct sum of 2 x 2 Mokhov's 
operator fT4l 

p _ ( ~^ yl dy ~ v l d x + v 2 d y + 2 u 2 \ 

\d x +v 2 dy-v 2 0 J 

and trivial 1x1 operator P = 0. 

Solution 2. The second solution is given by 

p = p(u>)W, 9 = 6i' = 6l‘=6| 1 = -6f = |, 


7,12 _ _ K 112 _ / 
0 2 K, °3 — q j 


where k ^ 0 is constant. This leads to 


P (y) = 


^{jP + KU 1 ) dy + 

(«-¥) ^ 


Kul+p'ul 


2 

KU?. 


\ 


KU y 

2 


(q - d y - KU 2 y + q’u 3 y 

0 0 

0 0 j 


The group (f50l) acts on this case as the previous one. Thus, we can reduce p and q to zero, and k to 
—2, obtaining 


(— 2 v 3 d y — v 1 v 2 d y + 2 v 2 v 3N 


P (v) = 


\ 


Solution 3. In the case where k = 0, the solution is given by 


11 p' 712 r 1 . 791 / 799 v' 


p = p(u 3 ), q = q(u 3 ), r = r{u 3 ), 6^ = £3 2 = v(u 3 ), bf = q' - v, bf = ^ 
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This leads to 


P (y) 


( p d y + 1 —^ L q d y + vu 3 0 
qdy + (q' - v)u 3 r dy + r -^ L 0 
0 0 0 


Here, the group d50l l acts on the objects as 

P -t P{p 3 )e~ 2 ' 4 ’, q -t q(<p 3 ), r ->• r(<p 3 )e 2 ^, v ->• (<p 3 )V(v? 3 ). 

Here we have four arbitrary functions p, q, r, v and only one function ip 3 and one constant i/j acting 
on them. Thus, even if we could consider several cases (where some functions are constant or 
zero), it does not simplify the classification. However, let us make a choice: if v is non-zero, it can 
be always reduced to 1. Thus the operator leads to 

/ 


P (y) = 


. p U . 

p dy H - 7f- q dy + eu : 


3 


o\ 


q d y + (q' - e)u 3 r d y + ^- 0 

0 0 ' 0 j 


(51) 


with e equal either to 0 or 1. 

Let us deal with d48t . If r = 0, solutions of conditions © lead to ( f5Tb (replacing r with 0). 
Otherwise, we have 


p = p(u 3 ), q = (u 3 ), r = r(u 3 ), b 1 3 1 = ^, b 12 = q ', b \ 3 = r', 


which leads to 


P (y) = 


v u y 


pdy+ 2 
qdy 


qd y + q'u 3 r d y + r'u 3 y ' 
0 0 


V 

In this case, © transform p, q, r as 


. - 2 , 3 , 3 , 2(3)V(3)\ , , . 

P e )-py- J > 9 “»• 9{<P 


0 


3, (P 2 )'r(p 3 ) 

)-—, r 


r(p 3 ) 

(<p 3 y 1 


,2b 


(v 3 y ’ 

Thus, since r ^ 0, we can always reduce r to 1 and p and q to 0, obtaining 

/ 0 0 dy 

P (y) = 0 0 0 

ydy 0 0 

Case iTTlL . Here the first structure is <EJ 2 , and the group of transformations which preserve its 
form is given by 

u 1 = e^V + <pV), u 2 = e-^V, u 3 = <p 3 (v 3 ). (52) 

In the case where the second metric is of the form ll47t . conditions |5| lead to two structures, 
given respectively by 


li _ P 


p = p(u), q = q(u), r = 0, b 3 = 




(53) 


and 


p = p(vr)u 1 + p(u 3 ), q = k, — 


Hd, r = 0 , = = jii = pV+£ U2 

2 1 ^ 2 ’3 2 


= -P, 


a- = 


i.31 _ E13 _ 

J °3 — °3 — 2 ’ 

u z 


(54) 
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where k is constant. 

In the case where the second structure is of the form given by d48l >. conditions <0 lead again to 
two structures, where the first is the same as (l53l> , and the second is given by 


p =p(u 3 )u 1 + p{u 3 ), q = K-- 1-, r = r(u 3 ), bf 1 = If = 


bl 1 = 


p'u 1 


h 12 

°2 


= ~P, 


bl 2 = - 


/ 2 
pu z 


hu — r _— 

°3 ~ °3 ~ r o’ 

U z 


— -h 6i — — 
°2 — °2 ~ o ’ 


which corresponds to J54t replacing r, ^2 3 ^3 3 >^2 1 ^3 1 With 


r = r(u 3 ), bf = -bf = ^, bf = -bf=r-- 2 . 

u z u z 

Finally, considering d49h conditions © imply 

P = 0, q = q(u 3 ), r = r(u 3 ), b { 3 = ~bf = b 31 = -b? = ±, bf=q', bf = r'. 

Summarising, this case leads to three different structures, namely 



(p 

d 1 p ’ u * 

n r\ 



( 0 

n r\ 

ru l~i u l \ 


a y ' 2 

q Uy 

U 2 


q (ly 


^0 

II 

q 

d y + q'u 3 

0 

0 

> P U) — 

q dy + q’u 3 

0 

r d y + r'u 3 


\ 

U 2 

0 

0 ) 


qu%-Tul 

\ U 2 

rd v 

0 / 


(55) 


/ / 1 I ~n 1 , pidy + {p'u 1 +p)u 3 ( pu 2 \ 1 2 . , ru 2 +(r’u 2 -K,)u 3 \ 

(l™ + p) dy H- ! -2- y - 2 J d V-P u y - 2 ~ r d v H --^-“ 


P U) = 


(»-¥) 




V 


T KU'Z—rU,. 

r d v + * v 


0 0 / 

(56) 

These are the more general solutions assuming the first operator given by (TTTll n. In these cases, a 
transformation of the form ll52t allows to simplify these structures, but we will get more cases. Let 
us discuss each operator in detail. 

Let us consider an operator of the form d55t 1 . Under change of coordinates of the form (l52l) . p 
and q transform as 

P -t e^p(<p 3 ), q -> q(<p 3 ). 

Thus, if p vanishes and q = £ is constant, we get g = fg, b = £&, that is, a trivial operator. If p 
vanishes, but q is arbitrary, it can be easily reduce to v 3 . Otherwise, if p ^ 0, it can be always 
reduced to 1, and, if q is not constant, the freedom in u 3 = p 3 (v 3 ) allows to reduce q to v 3 . Therefore, 
in this case we get the following non-trivial canonical forms 


P (y) 



V 3 dy 

3 3 \ 

_ v v v \ 

V 2 


f d y 

K dy 

_ Kv v\ 
v ' 2 

v 3 dy + vl 

0 

0 

’ P (y) = 

K dy 

0 

0 

\ V 2 

0 

0 J 

KV l 

\ V 2 

0 

0 J 


where n is constant and e can be either 0 or 1. 


When the operator takes the form (l55lh . if r vanishes, it reduces to the first operator of the 
previous case with e = 0. In general, a change of coordinates given by l l52t with if) = const, 
transforms r and q as 


(V? 3 ) 

(p 3 )' ’ 


q -> q(<p 3 ) 


( i p 3 y 
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Notice that here we have to impose the constraint ip = const, otherwise we would have an extra 
function in the metric written in the new coordinates. Thus, for r ^ 0, choosing ip 1 , <p 3 such that 
(p 3 )' = r(<p 3 ), (p 1 )' = <?(p 3 ), r can be brought to 1 and q to 0 , obtaining 


P (y) ~ 


0 Si 

0 d 


dy 0 


Finally, we have to look at the case (l56b . In the general case, a change of coordinates of the form 
transforms the functions p, p and r as 

r(p 3 )e-^ 


V -»• ( VW 3 ) - ) e p ^ (p(p 3 ) +?W - ) ) e 2 ^, 


(p 3 )' 


Thus, if r 7 ^ 0, we can reduce p and p to zero and r to 1. Otherwise, if r = 0 and p / 0, we can 
brought p to 1 and p to zero. If both r and p are equal to zero, then p can be reduced to 1. Finally, if 
r = p = p = 0, we have g = Kg, b = nb. Therefore, the canonical form of d56l can be summarised as 
follow 


( 0 

K dy 

K dy 

0 

J 1 v l~ Kv3 y\ 

a y H-^ — 

0 

0P 

se 

II 

( dy 

K dy 

K dy 

0 

_ KV v\ 

V 2 

0 

3 2 

J , Kv v~ v v 
\a y H-ps — 

0 

0 J 

\ ~v r 

0 

0 ) 


P (y) — 


t 


V 1 dy + 





0 


0 / 


Remark. Allowing linear change of the independent variables x and y, one can easily see that it is 
always possible to reduce to zero the part of P (y) which is proportional to Pm. This means that we 
can set k equal to zero in the above canonical forms. 


Case ITTh . Here the first structure is given by {TT} ;i . In the first case, that is, when the second metric 
is given by d47b . conditions (O imply 


bf = -p- 


u 3 p' 



112 «V 

f 3 ~ 2 , 

II 

CO 

t-h CO 
1-0 

II 

coco 

1-0 

i? = u* v+ ^\ 

bf = ~bf = 

u 3 (2u 3 p — k) 


2 u 3 p — i 


Thus, the operator leads 

/ 


P (y) = 


Pdy+ ^ 


(k — u 3 p) d v — 

(2 U 3 p—K 


(2 p+u p )u 


(K - U 3 p) dy - U -^ 

{U 3 f P dy ‘ (2 “ 3P+( “ 3 ) 2P> “ 


*y 1 2 

u 3 (2u 3 p—K,)u 3 


_ ( 2u 3 p-K,)ul \ 
u^u 1 —u 2 
u 3 (2u 3 p— K,)u 3 
u^u 1 —u 2 

0 


where p = p(u 3 ) and ft is constant. First of all, by linear change of x and y, we can set ft to 0. A 
change of coordinates which preserves P( x )/ namely 


u 


1 



\/v(v 3 )' 



cy/ip(v 3 ), u 3 


<p{v 3 ), 
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where c = const, transforms p into ' pp ^ . This means that, if p f -4 for £ constant, we can always 
reduce p to 1 , obtaining 


P (v) = 


~ V 3 dy - v 3 


— V 3 dy 

(v 3 ) 2 dy + V 3 Vy 

2 (V 3 ) 2 Vy 


2v 3 v 3 
v^v 1 —v'“ 

U^'/jl -11^ 


Otherwise, setting p = dL, we get 



/ Oly 

u 3 

(u 3 ) 2 

-Pd iiti 

^ 2 u 3 

2 «< \ 
u^u 1 —u 2 

P (v) = 

—l; dy 

_ 

2 u 6 

£u 3 d y —1— — 

2 £u 3 u 3 y 

—5-T- T 

u°u L —u* 


1 2Hu 

\ u^u 1 - 

y 

-u 2 

2£ n 3 u 3 v 
u^u 1 — u 2 

0 / 


In the cases where the second metric is given by ll48l l or by d49l) . conditions 101 imply 


p = r = 0 , q = k, bl 3 = — 6? 1 = —5 - w, b% 3 = —bo 2 — -=- 

u L u 6 — ir u L i 

where k is constant. This leads to a trivial operator, since g = ng and b = nb. 
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